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Abstract 

The Laplace operator acting on antisymmetric tensor fields in a D-dimensional 
Euclidean ball is studied. Gauge-invariant local boundary conditions (absolute and 
relative ones, in the language of Gilkey) are considered. The eigenfuctions of the 

•rH . 

^ . operator are found explicitly for all values of D. Using in a row a number of basic 

techniques, as Mellin transforms, deformation and shifting of the complex integration 
contour, and pole compensation, the zeta function of the operator is obtained. From 
its expression, in particular, ("(0) and C'(0) are evaluated exactly. A table is given in 
the paper for D = 3,4, ...,8. The functional determinants and Casimir energies are 
obtained for D = 3,4, 6. 
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1 Introduction 



In this paper we obtain the zeta function of the Laplace operator acting on antisymmet- 
ric tensor fields defined in a D- dimensional ball with gauge-invariant boundary conditions. 
Mathematically this computation is quite an imposing challenge, as is proven by the num- 
ber of erroneous results reported in the literature on this and related computations (details 
will be given later). The physical motivations for such a study are to be found in quantum 
cosmology, where the ( function of the Laplacian describes the contribution of antisymmet- 
ric tensor fields and ghosts to the pre-factor of the wave function of the universe (see e.g. 
Ref. [T|]). An intriguing problem in this context is the non- compensation of the boundary 
contributions to the one-loop divergences between different members of the supergravity su- 
permultiplet @. Another motivation is to provide the numerical material needed to extend 
previous analysis of the heat kernel asymptotics to the case of mixed boundary conditions. 

There are two admisible sets of gauge-invariant local boundary conditions — which have 
been called by Gilkey, respectively, absolute and relative boundary conditions |4|]. These 
sets are dual to each other and are becoming highly interesting in connection with recent 



developments in string theory. Hence one can study p-forms with p < 



D+l 



for both 



2 

types of boundary conditions. Due to duality, the determinant of the Laplacian for p-forms 
with absolute boundary conditions is the same as for (D — p)-forms with relative boundary 
conditions. Furthermore, owing to gauge-invariance we can restrict ourselves to transversal 
p-forms. The complete result will just be a sum of the ones for transversal p- and (p — In- 
forms, provided zero modes are properly taken into account. To find the spectrum of the 
Laplace operator on transversal p-forms we use the powerful method proposed in Refs. H 
They involve integral representations of the spectral sum, Mellin transformations, non-trivial 
commutation of series and integrals and skillful analytic continuation of zeta functions on 
the complex plane. Here we will focus again on a class of situations for which the eigenvalues 
of the differential operator, A, are not known explicitly, but where nevertheless the exact 
calculation of In det A is possible. The method is applicable whenever an implicit equation 
satisfied by the eigenvalues is at hand and some asymptoticity properties of the equation are 
known too. 

In this paper, we will find explicit solutions of the D-dimensional transversality condition 
in terms of p- and (p — l)-forms obeying a (D — 1) -dimensional transversality condition. 
These forms will satisfy now pure Dirichlet or Robin boundary conditions, instead of mixed 
absolute or relative boundary conditions. This clever procedure will enable us to find exact 
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eigenf unctions and to express the eigenvalues of the Laplace operator in terms of Bessel func- 
tions and their derivatives. After that, we will be able to perform explicitly the evaluation of 
the zeta function at the origin, and of the determinant of the Laplacian as well. We will also 
calculate the Casimir energy. A table of results is given in the paper for D = 3,4, 5, 6, 7, 8, 
which cover the situations which appear in the usual supersymmetric theories. However, 
our expressions to come are actually valid (and can be used) for any dimension D and yield 
explicit, exact values in a reasonable amount of algebraic computation time. Usual methods 
for the acceleration of the series convergence improve performance considerably. 

In connection with previous results, we should point out that for some scattered values of 
D = 4 and p = 1,2, several first heat-kernel coefficients have been calculated in the papers 
IH"0- These results agree with the analytical formulas in || once the corrections that were 
found in Ref. f| are taken into account (see also Ref. f2|). For D — 4, p — 1, the one- 
loop effective action has been evaluated in || for a specific choice of gauge and boundary 
conditions. 

The paper is organized as follows. In Sect. 2 we use the Hodge-de Rham decomposition 
of p-forms in order to simplify the structure of the spectrum of the Laplacian operator in 
the D-dimensional ball and, subsequently, of the corresponding determinant. After writing 
the absolute and relative boundary conditions as Dirichlet and Robin ones, a convenient 
analytical continuation of the associated zeta function is performed in Sect. 3, in some 
detail, what leaves us in a position wherefrom we can find the heat-kernel coefficients, the 
determinant, Casimir energies, and so on. The calculation of the zeta function at the origin is 
undertaken in Sect. 4, that of the determinant in Sect. 5, and the Casimir energy is obtained 
in Sect. 6. Finally, in an Appendix we give an exhaustive list of all the determinants obtained, 
both for the case of absolute and relative boundary conditions. 

2 Spectrum of the Laplace operator in a ball 

Consider the D = d + 1 dimensional unit disk with the metric 

ds 2 = dr 2 + r 2 dtt 2 , < r < 1, (1) 

where dfl 2 is the metric on the unit sphere S d . Throughout this paper we shall use the 
notations {x^}= {xq, Xi}, x° — r, — 0, 1, d. The (d + l)-dimensional Laplace operator 



3 



A = V M V M acting on a p-form, B, can be written as 



"V'B,,...,^ (2) 
(AB),,.., = (fl8 + + ^ + ""A) *m, + 

+ — ^ Vi a .Bj 1 ...j a _ 1 oi a+1 ...i p , (3) 

r 0=1 

where and "A are the covariant derivative and Laplace operator corresponding to 

the (i-dimensional metric g^. 

Any p-form B p admits the Hodge-de Rham decomposition: 

B p ^ B P± + dB ( P -i)± j ( 4 ) 

where B p± denotes a transversal p-form. The decomposition (^) commutes with the Laplace 
operator. Thus, in order to define the spectrum of the Laplacian on the space of all anti- 
symmetric forms, it is enough to study the case of transversal forms only. 

There are two sets of local boundary conditions consistent with the decomposition (|]). 
They are the so-called absolute and relative boundary conditions ||. In the coordinate 
system ([]]) the absolute boundary conditions read 

<9o-Bu,...,ijJ,9Af = 0, -Bo,ii,...,i p _il<9Af = 0, (5) 

while the relative boundary conditions have the form 

Bi lr ..j p \dM = 0, ^d H ^ J B 0jily „ t i p _ 1 \g M = 0. (6) 

Consider the (d + l)-dimensional transversality condition 

VB^ = 0. (7) 

On a disk it can be written as 

(VB) h ... ip _ 20 = V>V i B iil ... ip _ 2 o = 0, (8) 
(V5),,. Vl = (d + f—)B oil .. Ap _ 1 + ^V^,.,^ = 0. (9) 

According to the general method developed in the papers H |), the solutions of the equations 
(|S]) and (|5p can be expressed in terms of (i-dimensional transversal forms: 

B p± = B pT + B p± (^ p ~ l)T ), (10) 



where such (i-dimensional transversal forms satisfy the equations: 

m vm£ Vl = o, a£ ip , = o. (ii) 

Here A T denotes either B T or ip T . The second term in ( |I0|) has the form: 

BL,a* T ) - (- w A + (p " 1)d ; (p " 1)2 )^u- 

BU(V) = (4 + ^)r%*L^| (12) 

P 

^[ii'0i2— <p] := — 1) + Vj n '0i 1 ...j n _ 1 i n+1 ...j p 
n=l 



One can prove that the Laplace operator (0), ([3]), commutes with the decomposition (10): 

AB P± = AB pT + B p± (A^ p -^ T ). (13) 

The determinant of the Laplace operator on the space of (d + l)-dimensional transversal 
p forms can be represented as a product of two determinants, taken over d- dimensional 
transversal p- and (p — l)-forms: 

det(-A) p± = det(-A) pT x det(-A) (p _ 1)T . (14) 

Moreover, the fields B T and tp T satisfy pure boundary conditions. The boundary conditions 
for B T are defined by the first equations in (|5]) and @, for absolute and relative boundary 
conditions, respectively. For absolute boundary conditions on the field B p± , the form ^ T ( p_1 ) 
satisfies Dirichlet boundary conditions 

V\dM = 0. (15) 
For relative boundary conditions we have, for the (p — l)-form ip T , 

(flb + d ' 2P r +l ) V\8M = 0, (16) 

that is, Robin boundary conditions. We thus see that the initial eigenvalue problem for the 
(d + l)-dimensional transversal p-forms with mixed boundary conditions is reduced to two 
eigenvalue problems, for (i-dimensional transversal p- and (p — l)-forms with pure boundary 
conditions (Dirichlet and Robin). 

In the particular cases when p — 1,2, the boundary conditions (|i~5f ) and (^) agree with 
the corresponding expressions of papers |5|, 0. Note that the I = scalar mode generates 
a zero mode of the mapping ip — > B 1X . Hence this mode should be excluded when one 
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considers the path integral over transversal 1-forms and from the second determinant on the 
r.h.s. of Eq. 

Let us introduce the set of d- dimensional spherical harmonics, Y^ p i p (xj), corresponding 
to transversal p-forms on S d . They are eigenmodes of the <i-dimensional Laplacian. The 
associated eigenvalues and degeneracies Df are found to be |12], [13 1 



= ^[-Kl + d-l)+ P ]Y^ ip (x j; , 



gz + d-i) (i + d-iy. 
1 p\(d-p-i)\(i-i)\(i+p)(i + d-p-iy 1 ' 

We can represent the eigenf unctions of the complete (D = d + l)-dimensional Laplace 
operator as a Fourier series in the harmonics (TiT 



Bt..^^) = Y. Y Sl^)f^)- ( 18 ) 

(0 

Here we need to sum over (/), which means summation over the index I from 1 to oo and 
over another index, from 1 to Df, which describes the different harmonics with degenerate 
eigenvalues of the Laplacian ^A. This last summation is not shown explicitly. 
We can now substitute the decomposition ( |18D in the eigenvalue equation 

^ A - C Vl,£t2,---,A«p fll,(ij,-.fp \ LV J 

for the (d + l)-dimensional Laplacian (Q), (|3]). Let us recall the fact that for the fields B pT 
the zero-th components vanish identically: d _ 1 = 0. The equation for the components 
(|2|) reduces to the trivial identity = 0. The other components (0) lead to an equation of 
Bessel type for f^ l \r). After a rather lengthy algebra, one finds that the eigenfunctions of 
the Laplace operator (§) have the following form: 

rN^j^^^), (20) 

where J n denote Bessel functions. The eigenvalues \j are defined by boundary conditions 
and their degeneracies Df are given by (|T7|) . 

^From the preceding expressions, we are able to evaluate the determinant of the Laplace 
operator on the space of transversal p-forms. We obtain 

deV(-A) = nA? D ' 114 k , (21) 

1=1 k=l 

where for absolute boundary conditions the eigenvalues A and k are defined by @ and (|l5l), 
namely 

aor (i-d)/ 2+PJ(d _ i)/2+i(v) | aM = 0) J {d _ 1)/2+k (K k r)\ dM = 0. (22) 



For relative boundary conditions we have, from (^) and fll6|), 

J( d -i)/2+i{\ir)\ dM = 0, (do + d - 2p + i) r (i-rf)/2+P-i J (d _ 1)/2+fc ( Kfe r)| SM = 0. (23) 

One can easily check that the eigenfunctions defined in this section satisfy all necessary 
orthogonality properties. 

3 Analytical continuation of the zeta function 

Both the absolute and the relative boundary conditions can be written as Dirichlet and 
Robin boundary conditions, in the form 

^(V)U = 0, (24) 

u(d,p)J u (\tr) + \Mkr)\ dM = 0. (25) 

For the absolute boundary conditions we have u(d,p) = (1 — d)/2 +p, while for the relative 
boundary condition, u(d,p) = (1 + d)/2 — p. The boundary DM is here described by r = a. 
The zeta function is 

oo 

C d P A s ) = T, D fW 2s - (26) 

1=1 

The decomposition described in the last section will, at the level of zeta functions, manifest 
itself as a sum of the different zeta functions belonging to each term of the decomposition. 
Thus we can write 

C d p ±(s) = C d pT (s) + Ci > _ 1)T (s). (27) 

The zeta function is in general convergent for s > ^±1 only, but it can be analytically 
continued in the complex plane to all values of s, in particular to the vicnity of s = 0. 
Several authors have considered zeta functions corresponding to operators which eigenvalues 
are not given explicitly. In special, they have investigated in detail the case when they are 
given under the form of roots of equations involving Bessel functions (see, for example [fTofl). 
In the Refs. [|K], |11|] it has been shown explicitly how this analytical continuation can be 



carried out for zeta functions of this kind, and we will follow them closely. The reader may 
resort to those papers for all particularities skipped from the present calculation. 

Writing the boundary conditions (p4] ) and (p5"l) symbolically as ^(A/r)| QM = 0, the first 
idea is to express the zeta function as a contour integral along a path 7 enclosing all positive 
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solutions of the boundary condition equation, namely 



= E A(P, d) jT ^ 2 + ™ 2 r- In ^(o*). 



(28) 



Here we have introduced the constant m in order to simplify the treatement of the problem. 
It is, however, not essential to obtain the final result and we shall let this constant vanish later 



in the calculation. To be able to use most of the techniques developed in the papers QTOj, [IT 
we expand the degeneracy as 



d-l 



Di(p,d) = 5> i (d,p) [1 + 



d-l\ j 



(29) 



3=0 



The zeta function reads then 

d-l 



d-l oo / j _ -i\J r jj, a 

&00 = E^,p)E M + V J^ 2 + m 2 r^M l+ ^(ak) 



j=0 
d-l 



1=1 

oo 



E^(ME(' + ^PJ / 7 £ (fc2 + m2) ~^ lnvI/ ^ (afc) - (30) 

Now we are in the position of doing the analytic continuation. This involves substracting 
and adding the leading asymptotic terms of the uniform expansion of the Bessel function 



I u {k) and its derivative. For v — > oo and z = kjv being fixed, they are fUfl : 



and 



1 



v / 2^7(l + ^) 1 / 4 
1 e^(l + ^ 2 ) 1 / 4 



1 + E 



itfc(t) 



k=l 



i + E 



fc=l 



(31) 



(32) 



V27TZ/ -2 

respectively. Here Uk and are functions obtained in a recursive way in Ul4 |, while t = 
1 / y/1 + z 2 and // = vl + + ln[z/(l + \/l + £ 2 )]- Furthermore, we define the coefficients 
and M„(t) by 



In 



1 + E 



u k (t) 



k=l 



and 



In 



i + E 



Vk{t) u(p, d) 

f„ "T" 



Mi + E 



<E 

n=l 
Mfc(t)' 



V" 



k=l 



k=l 



E 

n=l 



M n (p, d)(t) 



V" 



(33) 



(34) 



Then, by adding and subtracting the first N terms of these last two expansions, we can write 
the zeta function for Dirichlet boundary conditions as 

N 



&00 = E A(s) + Z N (s) 



(35) 



i=-l 
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where, with v = I + and m = 0, we have 



o 2s T(q — I) d-l j, I 

^-i(') = , f4 . 2 A E e^lOU* - 1 - j, ^), 



40Fr(s + l) 



,2s d-1 



Ms) 



■°— ]T ] (d,p) (h (2s - j, — ) , 
4 i=o z 



,2s d-1 



2r(a) 



fc=0 



r(i + fc + 



and 



(36) 

(37) 
(38) 



Zm(s] 



2sa 



2.s 



Sill 7TS 



7T 



j -2s 



3=0 



Z=0 



x / cfez 



-2s-l 



In Iv(vz) — In 



n=l 



(39) 



.v / 2^(l + ^ 2 ) 1/4 

The coefficients x k ,i in equation fl38"|) are obtained from the polynomial expansion of Di(t): 



k=0 

Similarly, denoting by the coefficients in the expansion of Mi(t), 



2 1 



Mi(p,d)(t) = J2zi, k (p,d)t 

k=0 



: +k 



we can write the zeta function for Robin boundary conditions as 

N 



i=-l 



where 



A* l {s) = A_ 1 {s), 
A*(s) = -Ms), 



A*(s) 



•^-y £ e j (d,p)Q H {2s + 1-3, — ^— ) 2^ z kti (p, d)- 



k=0 



r(i + *±*) 



and 



(40) 



(41) 



(42) 

(43) 
(44) 

(45) 



2sa 2s b -^llL V eAd,p) V u j ~ 2s / dzz^H In [u(p, d)Uvz) + zvI'Jvz)] 

71 j=0 1=0 JO [ 



-In 



v 



e^(l + z 2 ) 1/4 



_ £ M n (t) 



n=l 



V" 



(46) 



It can be shown that both zeta functions, ( |35"D and fl42|), are well defined for d ~ l ~ N < s. 
Both Zjy(s) and Z§(s) are analytic here, so that all the poles are contained in the A's. The 
analytical continuation can therefore reach the desired range of s, by just changing the value 
of N. We are thus in a position where we can find the heat kernel coefficients, the zeta 
function determinant and also the Casimir energy. 

In the following we shall work in the unit sphere. Since the only change we have to do 
on this zeta function in order to include an arbitrary radius is to multiply by a 2s , the results 
that we obtain for the unit sphere can be easily converted into those for the general case. 



4 Calculation of the zeta function at s = 

^From the zeta functions we have defined above, we can now calculate the heat kernel 
coefficients, using the relations that exist betwen them. ^From a physical (and maybe 
also from a mathematical) point of view the most interesting coefficient is the one that 
corresponds to C(0): 

C(0) = -in, (47) 
where the numerator Bd+i is the corresponding coefficient that comes from the short time 

2 

expansion of the integrated heat kernel: 

oo 

K{t) ~ (4vrt)^ £ Brutf. (48) 
m=0 

Note that Zd(0) = and Z^(0) = 0, since the sum over I and integral over z is convergent 
here for N = d. Therefore, we need only consider Ai(0), i = —1, 0, . . . d. Using the expansions 

(-1)" 

T{-n + e) ~ L -r- + .. (49) 
n\e 

and 

C(l + £) z/) --?--*(*/), (50) 

we find the expressions 

A-i(0) = ^£ ei (d,p)C*(-l -j,^-), (51) 

4>(0) = -j'zZ^MpXhH,^), (52) 

j=0 
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v\ d 


7 


6 


5 


4 


3 


2 


A 


3559 
9072 












3 


36583 
45360 


2929 
4608 


358 
945 








2 


20467 
25200 


1624993 
1935360 


1531 
1890 


81 
128 


7 

20 




1 


185449 
226800 


785567 
967680 


6199 
7560 


2429 
2880 


49 
60 


5 
8 





3629089 
3628800 


1934993 
1935360 


6379 
7560 


9647 
11520 


151 
180 


41 
48 



Table 1: £(0) for absolute boundary conditions, transversal p-forms. 



p\d 


7 


6 


5 


4 


3 


2 


3 


3559 
9072 


2929 
4608 










2 


4283 
25200 


416417 
1935360 


358 
945 


81 
128 






1 


33521 
226800 


143263 
967680 


1109 
7560 


541 
2880 


7 
20 


5 
8 





289 


367 


1 


17 


1 


1 


3628800 


1935360 


1512 


11520 


180 


48 



Table 2: ((0) for relative boundary conditions, transversal p-forms. 



Mo) = --e^ XX* r(i + fc + i) ' (53) 

and 

A R (o) = -KiE^(p^)^^tw- ( 54 ) 

* k=0 1 \ l ~r 2 ) 

The numerical values obtained from these expressions are given in Tables [I] and j|. Only 
those values that are independent have been given (the rest are got using duality). For 
absolute boundary conditions and p — 1, the scalar field B 0± has a zero mode satisfying the 
boundary condition doB 0± = 0. But this mode does not contribute in the Hodge-de Rham 
decomposition. The scalar field is treated by construction of the zeta function for Neumann 
boundary condition. From the definition of the integrated heat kernel 

K(t) = J2e- X "\ (55) 

n 

we see that omision the zero mode corresponds to subtracting 1 from this sum. Since the 
relationship (|47|) is still valid — also when the zero mode is projected out — we conclude that 



11 



in order to get ((0) without zero mode we need to consider Q nc i{0) — 1. In this way we have 
obtained the values also for p = 0, by extracting the values of Cincz(O) from |l(J. For absolute 
boundary conditions and d = 3, all values are in agreement with those calculated in 0. For 
relative boundary conditions and d — 3, p — 1, we have found the value given in H. 



5 Calculation of the determinants 



We will employ the zeta function definition of the determinant of an operator, A, namely 

d 

- C-'(.s) 

s=0 



lndet(A) = -^-C A (^ 



(56) 



The determinant of our Laplacian is thus 

-lndet^(A) = C^(0) + C ( ti)T(0) 



(57) 



The determinant for the sphere of radius a is obtained by adding the terms coming over 
from the derivation of a 2s : 

- lndet^(A)(a) = 21na [C d pT (0) + C^ 1)r (0)] + $.(0) + $-i)t(0). (58) 
The differentiation of the A's can be done quite easily. Following the same steps as in 



Ref. IT], one obtains the formulas 

d-i 



z' d {o,x) = J2 e Mp){ / dtt \ 



^) d? ( r 1 



e~* dP \ e* - 1 



* D n (l)T(x + n-j) , . d + 1 

+ E — T\ ( H [x + n-j, ) 

n=i+l 



r(n-j) 



'-iy(j + i)\( H ( x -j-i, 



d+l. 



T(x-j-l) 



and 



< j+l U (p,dy + i d+i. 



Z' d (0,x) + Y / e j (d,P)\(-^ 

3=0 I 

+ e ^^(-ir +1 c^ — 



d + 1. 



n=j+2 



n 



,d + 1 



+ a:) 



-f-lV'j [ U dxx^HnT 
Jo 

+(-l) j u(p, d) \nT( d ^ + u(p, d))\. 



(59) 



(60) 
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The parameter x is introduced in order to allow for the individual terms to be finite. In the 
final answer this parameter will disappear. The determinants obtained in this way are listed 
in the Appendix. We have also included the determinants for transversal p = forms given 
in Ref. JTTJ. For Neumann boundary conditions the zero mode must be treated specially, 
what yelds an answer we can use directly later on. 



6 The Casimir energy 

As is well known, the Casimir energy (or vacuum energy) density can be written as a (usually 
formal) sum over the eigenvalues of the energy equation, that is | J2k ■ The energy density 
difference gives rise to the Casimir force. However, this sum is usually divergent and has to be 
regularized. A very simple and elegant way of performing the regularization is to use the zeta 
function method (see |15|, |TjJ for very extense and updated expositions of this procedure). 
But, sometimes, it happens that even after analytical continuation the zeta function at the 
desired value still diverges. The normal procedure consists then in resorting to the principal 
part prescription flTT| , |T8|j . In Ref. [[18] the physical meaning of this prescription has been 
investigated in depth. A finite part of the vacuum energy is found by separating off the pole. 
Obviosly, from our zeta function, the vacuum energy is obtained by computing its value at 
s = —1/2. Writing the zeta function around s = —1/2 



(61) 



we see that the vacuum energy is given by 



Ec = (62) 

At s = —1/2 one observes that the poles come from the gamma function of A_x, from the 
gamma function of Aj, for b = and i — 1, and from the zeta function of A$, when i = m + 2. 
We perform a Laurent expansion around these poles and isolate the corresponding finite 
parts. The rest of the functions will only contribute to the finite part. The values Zd+i(—h) 
and Zj +1 (— ~) have to be computed numerically. These contributions are generally quite 
small compared with the finite part which comes from Y^f=-\^i- By adjusting N, the 
values of Z d+l {—^) can be further reduced, allowing one to obtain the same accuracy with 
much less effort. For some values of d and p the sum over I converges very slowly. Use 
of Richardson extrapolation leads to a dramatical improvement of the convergence speed. 
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d 


P 


absolute boundary conditions 


relative boundary conditions 


2 





0.008891 + Vi n i a) 


0.02806 + 2 4° (a) 




1 


0.1678 + 16 i" (a) 

315 7T 


0.1678 + 16 '" (a) 

315 7T 


3 





-0.001793 


-0.03537 




1 


n Q/ifio i 1339 ln ( a ) 
u.oiuz -t- 3276g 


0.04881 m ll^ a) 

65536 




2 


0.04881 16 ^' n fi (a) 

65536 




4 





0.000945 3 /'^ a) 

45045 7r 


0.03054 + 23 Z^ {a) 

15015 7t 




1 


„ . C r,r, , 11048 ln(a) 

0.4677+ 45045 ^ 


0.01881 + 66 i 2 J° (a) 

45U45 7r 




2 


0.1749 21 4 1 ' l(a) 

3465 7r 


0.1749 212 J' l{a) 

3465 7r 


5 





0.0002050 


n (MQ11 3118613 ln(a) 
U.UZ01Z 50331648 




1 


n^n | 871339 ln(a) 
u,di ^ J ~ 8388608 


n mrm 1052991 ln ( a ) 

U.UZU^I 16777216 




2 


no^rn 1063379 ln(o) 
yj.o^tov 25165824 


104 8576 




3 


n 0^70 i 31697 ln ( a ) 

U.U00M+ 1048576 





Table 3: for (d+ 1) -dimensional transversal p-forms. 



d/p 


5 


4 


3 


2 


1 


2 








0.1959 


0.1767 


3 






-0.08417 


0.2974 


0.3444 


4 




0.04935 


-0.1561 


0.2928 


0.4668 


5 


-0.04340 


0.03546 


-0.2902 


0.1790 


0.5251 



Table 4: (j) for p-forms on the unit sphere with absolute boundary conditions. 
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This extrapolation is a general procedure of numerical analysis. It is here valid because the 
partial sum has the following asymptotic behavior: 

]>> J+1 / {...}^Qo + Qin- 1 + Q 2 n- 2 + Q 3 n- 3 + ... : n^oo. (63) 
1=0 Jo 

The finite contributions for the (d + l)-transversal forms are listed in Table 3. We have 
included the scalar field, p — 0. When the argument of the zeta function is negative, the 
constant term for Neumann boundary conditions does not contribute. Special care of this 
term need therefore not be taken here. The coefficients belonging to ln(a) equal the heat- 

kernel coefficient — . Table 4 gives the vacuum energy for the unit sphere for all 

p-forms, for absolute boundary condition. For p = 1 we see that the energy increases with 
increasing d. For p = 2 there is a maximum at d — 3, while the energy for p = 3 and p = 4 
decreases with d. For constant d there are actually less systematic trends. 
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A The zeta function determinants 



A.l Absolute boundary conditions 

In this case, we have obtained 

, , c / a \ 55120073 4 r 1 , , T . , 1 r 1 , , T . 

-hided, A = / dyy 2 ty(3 + y) + - dyy 4 ^(3 + y) 

3±v ; 64864800 3 Jo V ; 3 Jo V ; 

29 1n2 215 1n3 + ln4 + CR(-5) Cr'(-4) 



3780 2 8 6 12 

_ Cr < ( _3) + z^ta + iMzii _ (Ao) 

_ lndetL( A) = -38814043 5179 1n2 _ 3h4 fa^-5) 

2±v ; 64864800 3780 8 6 

i Cr (-4) t u 7Cr / (-2) 4Cr / (-1) , 

- Cr (-3) ^— + + Cr (0) 
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-In detf ± (A) 



-lndeto ± (A) 



-In det2 ± (A) 

-In det? ± (A) 
-In deto ± (A) 

-In detij_(A) 
-In deto ± (A) 



75711793 8 



+ | / 1 dyylnT(3 + y)-^ f * dyy 3 lnT(3 + y) 
6 Jo o Jo 



64864800 

3 hi(|) _ 11869 In 2 _ Cr(~5) CrM) 

2 3780 12 8 

Cr'(-3) 5Cr'(-2) Cr'(-I) 



3 

7087979 
32432400 

2 



1181 In 2 + + Cfl(-5) + Cb'(-4) C*'(-2) Cfl(-l) 



3780 



60 



24 



4ndet^(A) = 



-lndet^(A) 



-lndet^A) 



~ jdy(y- 2) \nT(y) + l - J dy {y - 2) 3 \nV{y) 

o o 

f 2 dyyV^ + y)- f 2 dyy 3 + y) - 
Jo A Jo I 



3411 9 
2560 + 4 

15Cr'(-4) + 21Cr'(-2) 



24 



45 In 2 

64 



60 



32 

10633 9 



16 



J Q ~ 2 dy lnl^ + y) + 3 jT' dyy 2 Inl^ + y) + ln(|) 



11520 4 

| 581 ln2 5Cr(-4) 7Cr'(-3) | 5Cr'(-2) | 13Cr'(-1) 



2880 16 
19261 713 In 2 



24 



8 



24 



46080 

3/2 



720 



+ In 5 - Hdi ' A) - 7Cfi/( ~ 3) - C * ( ~ 2) + ("I) 



64 



48 



32 



48 



3/2 





+ C dyy 2 m{2 + y)- 
Jo 

+ Cr(0) 



o 

2081 
3360 
3Cr'(-1) 



In 2 ln4 Cr(-2) 

^o""~ + Cr( " 3) "^ — 



5989 83 1n2 3 In 4 , Cr'(~2) 3Cr / (-1) . 

H ~ 1" CR 1-3J H r ~ U (Oj 



10080 



60 



8 



2 
i 



493 61 In 2 C?(-3) Cr(-2) Cr(— 1) , / , 



4320 ' 90 3 2 6 

19 .,.3 3 . . 3 



_ + 2fA,l„r(- + ! ,) + i ln2--CK'(-2) 

71„2 . . SO^.Cs^l).,^^ 



7 

32 



1/2 



6 



+ ln3 - 



(64) 



A. 2 Relative boundary conditions 

In this case, the results are 

, , 5 /AX 55120073 4 /-i 2t/o , 1 r 1 , , T , n 
-lndetL(A) = P/)OP/)onn -- / dyy 2 V(3 + y) + - dyy 4 *(3 + y) 

3 jo 3 jo 



64864800 



16 



-lndetfj_(A) 



-lndeto ± (A) = 

-lndet^(A) = 

-lndet^(A) = 

-lndetJ ± (A) = 

-lndeti ± (A) = 

-lndeto ± (A) = 

-lndetij_(A) = 

-lndeto ± (A) = 



29 1n2 215 In 3 , In 4 , Cr'(~5) Cr'(~4) r l( 
+ — — H — Cr 



3780 2 
7Cr'(-2) , 4Cr'(-1) 



6 



12 



+ 



12 3 
-11173163 1 



Cr'(0) 



dyy 2 V(3 + y) + — cW^(3 + y)-- — ^ 
64864800 12 7o yy K y> 12 Jo yy K y> 2 



8341 In 2 , Cr'(-5) Cr'M) Cr'M) 



+ 



3780 12 
5Cr'(-2) Cr'(-I) 



4027 11 1 1 1 

'6486480 " 756 2 + 60 C -'- 5) " 24 C «<" 4 » + 24 C «<" 2 > " 60 ^ 



3411 9 



2560 



+ j jj dyy^^ + y) - jj dyy 3 q(^ + y) 



45 In 2 
64 



15Cr'(-4) + 21Cr'(-2) 



32 



16 



1 r^- 

1 / 2 



17021 1 /■§ , _ .5 

11520 + 12 Vo ^^"sio Vl/! 2 
5 Cr (-4) + 7Cr(-3) + 5 Cr'(-2) 13 Cr'(-I) 



2561 In 2 
2880 



24 



24 



16 

47 17 5 7 1 1 

— - + - ! — In 2 + — Cr(-4) + — Cr(-3) - — Cb(-2) - — C*(~l) 
9216 2880 64 Sin ' 48 ^ ; 32 SRV ; 48 RV ; 

2081 f 1 , 2 T , rt , ln2 ln4 , ., n , 
+ / dyy 2 V(2 + y)- — - — + Cr -3 
Jo 



3360 

Cr'(-2) 3Cr'(-1) 



20 



+ Cr'(0) 



_iZ 3 _ + l ln2 + I C '(- 3 ) 
30240 90 3 SRV ; 

19 
16 



+ 2 jT'dy lnr(^ + y) + ^ln2-^CR'(-2) 



■|-^n2-^(-2) + ^(-D 



(65) 
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